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1.0 INTRODUCTION 



1.1 THE PROBLEM UNDER CONSIDERATION 

It is the nature of submarine warfare to avoid detec- 
tion. Indeed, one of the major advantages of the submarine 
is surprise. The passive detection of low level acoustic 
signals (in a water environment) is a function of many 
factors. Among these are the power density of the signal 
at the detector and the period of time over which the signal 
is received. This paper seeks to examine the differences 
in suggested evasion tactics which result from a study of 
these two factors. 

The power density received by the detector may be 
represented as shown in Figure 1.1.1. 



Power 

density 

received 
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Due to the nature of the detection equipment, a memory 



or integration time is required. The integration time is 
shown as T in the figure. For an evader, which is producing 
the noise signal and wishes to minimize the probability of 
detection (all other factors being held constant) , the 
objective is to minimize the maximum area over time under 
the power density curve over the integration time. The 
memory requirement of the equipment is the same for each 
instant of time. While time advances this "window" of 
integration time, T, slides along with the power density 
curve. Expressed mathematically the objective is to: 



A solution in the general case is not proposed, only very 
large and very small values of T will be examined. That is, 
the extreme effects of integration time on the optimal 
evasive trajectory are studied. 

Typically for an evader at 60 miles range initially 
and a relative closing velocity of 20 knots an engagement 
with a searcher will be approximately three hours. 




T 



Min 

P 



P (u) du 



(i.i.D 



9 



Integration times of less than an hour are available with 
today's equipment. As an upper bound the case where the 
integration time is much greater than the engagement time 
may essentially be called infinite. The lower bound would 
be zero. The boundary cases are then "infinite" and "zero" 
integration time. The goal is to solve for the two optimal 
evasive trajectories and to compare them. 

The infinite integration time problem can be formulated 
to fit within the general framework of the control problem. 
This class of problems is an outgrowth of the calculus of 
variations. It is one objective of this paper to apply the 
maximum principle to the infinite integration time problem. 
As the integration time reduces to zero, the objective 
becomes to minimize the peak power. This problem is 
formulated in ordinary differential equations. 
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1.2 NOISE POWER RADIATION AND ATTENUATION 



The noise power radiation, for a line spectrum of a 
submarine has been modeled to have the mathematical form 



where; a, b are constant coefficients to be determined, 
and v is velocity, measured in knots. 

Data was obtained on the radiated noise levels of 
submarines as a function of velocity at a particular line 
frequency. A least squares curve fit was made to the data. 
The form of equation (1.2.1) was used. A correlation 

o 

coefficient, R , greater than .9 was obtained. This was 
considered sufficient justification for the use of equation 
(1.2.1) as a model. The power equation was then normalized 
to the form, 



The original data used for the development of the coeffi- 
cients is classified. Since the units in which power is 
is measured are unimportant as far as the analysis is con- 
cerned, and so that this paper could be unclassified, the 
coefficients were normalized. The final form of the power 
equation used was, 



P(v) = a + b V 4 



( 1 . 2 . 1 ) 



P(v) = 1.0 + b/a v 4 . 



( 1 . 2 . 2 ) 
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P(v) = 1.0 + 0.00005 v 4 . 



(1.2.3) 



A further assumption is that the submarine is an iso- 
tropic radiator. That is, the effects of lobe radiation 
are not considered. The attenuation of the noise power is 
also assumed to follow an inverse square law. This was 
strictly for convenience, since the solution can easily 
accomodate any attenuation function as long as it is a 
function of range. The simplifying assumption about iso- 
tropic radiation is more serious, and as yet, no manner of 
incorporation into the proposed solution has been discovered. 
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1.3 HYPOTHETICAL SUBMARINE ENCOUNTER 

Suppose submarine E, the evader, is assigned a patrol 
area in which she is motivated to remain, but may maneuver 
at any course and speed she chooses. Suppose further that 
submarine S is searching along a straight path with a 
constant search velocity V g for E. Assume also that E is 
directly on the search path of S as in Figure 1.3.1. 




Figure 1.3.1 



In order to attempt to avoid detection, E must detect 
S and then move away. The situation where E detects S first 
is only considered here. Once S detects E, the game is over. 
The hope is that E will know about the presence of S, and 
will be able to maneuver in such a manner that minimizes 
her probability of detection. The evader is not constrained 
in this problem by acceleration, and may therefore take on 
any initial velocity desired. The question is then, what 
is the optimal trajectory which minimizes the probability 
of detection? 
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The nature of the solution to this problem gives the 
answer to a set of possible initial evasion states. It 
is not necessary that the evader start from a point directly 
on the path of the searcher. Evasion may begin at any state 
position. 
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2.0 PROBLEM FORMULATION 



2.1 INFINITE INTEGRATION TIME FORMULATION 

As an extreme the integration time of the searcher* s 
equipment could be infinity. That is, for this situation 
the objective of the evader would be to minimize total noise 
energy received by the searcher over time. The total energy 
received by the searcher can be expressed as an integral: 



The power received by the searcher is a function of the 
evader* s velocity and distance. Assuming an inverse square 
law, 



i 



co 



total energy received=J = 



Power received (t) dt. 



o 



( 2 . 1 . 1 ) 



Power Received 




( 2 . 1 . 2 ) 



R 



The problem may now be expressed as, 




.oo 



(2.1.3) 



o 



where v and R are both functions of time. 
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This class of problems, which results from the calculus 
of variations, may be solved by one of two modern approaches 
to the control problem. (2) 

Older methods of static optimization involved the intro- 
duction of Lagrange multipliers, one for each constraint; 
defining a Lagrangian expression and finding a saddle point 
of this expression. 

One modem approach to the control problem is dynamic 
programing and the application of the principle of optimality, 
which states that: 

n An optimal policy has the property that, whatever 
the initial state and decision are, the remaining 
decisions must constitute an optimal policy with 
regard to the state resulting from the first 
decision." (3) 

The dynamic programming approach basically assumes an 
optimal performance function exists and defines a recurrance 
relationship. The partial derivative of this recurrance 
relationship is called Bellman's equation, the solution 
being the optimal performance function sought. 

Another modem approach is the maximum principle. (2) 

The maximum principle can be considered as an extension of 
the method of Lagrange multipliers to dynamic optimization 
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problems. The general control problem to be solved by the 
maximum principle is; 



Max J = I(x,u,t) dt + F(x-pt^) 




subject to: ^ = f(x,u,t) 



x < t o ) =x o 

x(t x ) =x 1 



( 2 . 1 . 4 ) 



[u(t )] £ U 



where I( ), F( ), and f( ) are continuously differentiable 
functions; tQ» x 0 are P arameters defining the terminal sur- 



the given control set U. 

In control theory there are state variables, control 
variables, and costate variables. In this problem the state 
variables, here represented by x are actually the x and y 
position coordinates, in relative space. The variables 
under "control" are here represented by u, but in the 
actual problem are velocity (v) and heading ( 0 ) . Velocity 
and course heading are both functions of time. 



face; and 




, the control trajectory must belong to 
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For each constraint a "costate" variable is introduced. 
This costate variable, A(t), is the dynamic equivalent of 
a Lagrange multiplier. This leads to the development of 
the Hamilitonian function which is defined to be; 

H(x,u,t) = I(x,u,t) + f(x,u,t). (2.1.5) 



The infinite integration time problem may now be expressed 
as ; 






Max J = 



- P(v) 

x 2 + y 2 



dt 



subject to; 



d x 
dt 



-V g + v sin 0 



= v 
dt 



cos 



0 



( 2 . 1 . 6 ) 



x(t 0 ) 

y(t 0 ) 

x(t 1 ) 

y(tp 



|> initial conditions 



> terminal conditions 



where, 0 is a function of time and is the course heading, 
and v, also a function of time, is the velocity. In order 
to achieve a practical mathematical solution, infinity is 
approximated with a large number. 
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The constraints dx/dt and dy/dt are obtained from the 
geometry of Figure 3.1.1. Note that the course heading, 0 , 
is measured from the vertical. 




Figure 2.1.1 - geometry relative to the searcher 
In this problem the correspondence of functions is; 



K 


, - P(v) 

y/1 + y2 


(2.1.7) 


F( 


) = 0 


(2.1.8) 


dx 

dt 


= - V s + v sin 0 


(2.1.9) 
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— ^ = v cos 0 



( 2 . 1 . 10 ) 



The solution to this problem is contained in appendix A. 

Just as the Lagrange multipliers of static problems 
yield information on the sensitivity of the solution, the 
costate variables of the maximum principle yield information 
on the sensitivity of the solution to variations in para- 
meters. 

Let J* be the optimal value of the objective function 

then 

■ a* <'0> ■ f 2 - 1 - 11 * 

That is, the sensitivity of the optimal objective function 
to changes in the initial state variables is given by the 
corresponding initial costate variables. If a costate 
variable is small then the relative sensitivity of the 
optimal objective function will be small for small changes 
of the initial state variable. 

In static problems the Lagrange multipliers have the 
interpretation of "shadow price" and a dynamic analogue 
exists . 
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2.1 ZERO INTEGRATION TIME FORMULATION AND SOLUTION 



As the integration time is reduced to zero, it can be 
argued that the objective of the evader is to minimize the 
peak power received by the searcher. The problem may be 
viewed as one in which there exists a threshold of power 
density level, which, if exceeded, would result in detection 
by the searcher. 

Assume there exists a set of state positions S where, 
for a given threshold of detection, the evader could not 
escape. Assume also there is another set S where the evader 
could always escape. The sets could take the form shown in 
Figure 2.2.1. 




There must exist some boundary between the two sets of state 
positions. It would be of interest to know the general 
shape of this boundary. 



21 



Suppose the evader is in some state as shown in Figure 
2.2.2. For a given power density threshold K there will be 
associated a corresponding velocity v k that satisfies, 



p (vk> 

x^ + y2 



K 



( 2 . 2 . 1 ) 



where v k is the maximum velocity the evader can have without 
exceeding the threshold of detection. Actually the evader 
is allowed to choose any velocity and heading within the 
circle of Figure 2.2.2. 




Figure 2.2.2 

The maximum velocity possible and the state position deter- 
mine the angle 0 k . That is, 

0 k = arc sin (v k /V g ) (2.2.2) 

By geometry 0 = 0 k . The question the evader must answer is, 
which controls should be chosen. 
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Suppose the evader is in set S, the set of states where 
he cannot escape detection no matter which controls are 
chosen. In Figure 2.2.3 it can be seen that the shape of 
the boundary must be such that its slope at the state 
position is greater than or equal to 0^. 




If the slope of the boundary were less than 0^ it would be 
possible for the evader to escape from the set S. This 
would be a contradiction. 

Suppose the evader is in the set of states where he 
can always escape, provided he chooses the right control. 
Since the evader is not allowed to choose a control which 
produces a 0 greater than 0^, the slope of the boundary 
must be less than or equal to 0^. Otherwise E would always 
move from the region of escape to non-escape as shown in 
Figure 2.2.4. 
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Figure 2.2.4 - boundary with slope greater than 0^ 

In order to Resolve these contradictions, the evader 
should choose the controls which produce a 0^ which is 
tangent to the boundary. This is shown in Figure 2.2.5. 




Figure 2.2.5 - boundary with slope equal to 0^. 

As shown by Figure 2.2.2 this 0^ is the result of choosing 
the maximum velocity v^. 

Since the choice of a velocity control will be always 
to maximize velocity for a given power level which means to 
always produce as much power as possible, the result is that 
the power level will remain constant at the detector. 
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For the zero integration time situation the objective 
is to maintain a constant power level received by the 
searcher. This may be expressed as 

P(vf) 

PFINAL = — 5 t : = K (2.2.3) 

*f + 7f 

where K is a constant, and the f subscript indicates the 
final value at the termination of the trajectory. 

The equations of motion may now be derived from the 
geometry as shown in Figure 2.2.6. 




Figure 2.2.6 



Velocity, v, and heading, 0 , are functions of time. The 
equations of motion are; 

= -V s +vsin0 (2.2.4) 

= v cos 0 (2.2.5) 
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Equations (2.2.4) and (2.2.5) may be combined to express 



£ = V C ° S 0 (2.2.6) 

dx - V g + v sin 0 

The problem now arises that the initial conditions 
which would allow a solution to the equations of motion are 
part of the results that are sought. One question that the 
evader would have, is what initial velocity and heading 
should be chosen. Further analysis of the nature of the 
trajectory should supply a solution. 

The following remarks refer to Figure 2.2.7. 




Suppose the evader is in state (xQ,y Q ) and has a given 
power level, K. This implies that v k (0) is greater than 
zero, where v k (0) is the velocity at state (x 0 ,y 0 ). There 
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will be associated with this power level, K, a trajectory, 
not necessarily optimal. 

Suppose the final velocity at x=0 is greater than zero. 
The game is considered over at x=0 since this would be the 
minimum range, and hence, the evader could simply stop all 
action and continue to improve his situation. It is con- 
tended that there exists a K’ less than K, since a v^i(O) 
could be chosen which is less than v^.(0) . This would imply 
that 



P (v k i(0)) ^ P (v k (0>) 

x 2 + 2 < x 2 + 2 • 

x 0 y 0 x 0 y 0 



(2.2.7) 



Since v^,(0) v^(0) , this implies that Initially, 

at least, the K 1 trajectory will lie below the K trajectory. 
Could K' cross K as indicated in the figure? If this were 
to occur, the 0^i would be greater than 0^, which would imply 
that K 1 is greater than K. This is a contradiction, and 
therefore, the trajectory for power level K 1 cannot cross 
that for K, but must lie below it. 

Since trajectory K' is below that for K, the terminating 
distance y^» must be less than y^. From equation (2.2.3) the 
conclusion is that v^ i must be less than v^. That is, the 
terminal velocity must decrease. The result then is to 
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choose a lower power level until the terminal velocity is 



zero. 



Could the terminal velocity go to zero before reaching 
the y axis? This is not a feasible result, since this is 
not the state position of minimum range. Due to the move- 
ment of the searcher the range would decrease and the power 
received would increase above the detection threshold. 

The object was to discover the initial conditions so 
the equations of motion could be solved. Instead the 
termination conditions are known. These will be used as 
initial conditions by reversing time. The system of equa- 
tions to be solved is now; 




( 2 . 2 . 8 ) 




( 2 . 2 . 9 ) 



and, 



x(t Q ) = 0 



y(t 0 ) ~ parameter 



v(t 0 ) = 0 
0(t o ) = 0 



are the initial conditions. 
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3.0 RESULTS 



The results of the analysis consist of giving advice 
on the choice of control variables to an evader. Both 
cases are presented; infinite integration time in graph 
(l.A.) and zero integration time in graph (l.B.) Advice 
is given in both cases only for the first quadrant, as this 
is the region of greatest concern. Four evasive trajectories 
were calculated and are shown on each graph. Those trajec- 
tories considered were for initial evasion ranges of 60, 

40, 30, and 20 miles. Trajectories were calculated for 
each of the four initial ranges both for infinite and zero 
integration time. The set of graphs labeled 2.-- are for 
a search velocity of 15 knots. Set 3.-- is a replication 
of set 2.-- but for a search velocity of 10 knots. 

For the 10 knot search velocity and infinite integra- 
tion time trajectories were calculated for initial evasion 
ranges of 30 and 20 miles . The 60 and 40 mile curves went 
off scale and were omitted. 

Control variable measurements were taken from the 
velocity and heading graphs and these control vectors were 
plotted on the relative movement plots in graphs l.A and 
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l.B. These vectors are the advice which would be given to 
an evader. It is not necessary that an evader start at a 
state position such that he is directly on the path of the 
searcher. He may start at any position and begin evasion 
with the control vector belonging to that state. If the 
searcher should change speed it would be necessary for the 
evader to consult a new set of control variables and continue 
the evasion. 

Graph l.C shows a comparison of the relative movement 
plots of infinite and zero integration times. Only the 60 
and 30 mile curves are presented for clarity. The control 
vectors for the two situations do not appear significantly 
different during the initial stages of the trajectories, 
and the result is similar trajectories. As the range de- 
creases the effects of integration time become more apparent. 
The graph shows that near the minimum ranges the effects of 
integration time are critical and hence it would be of 
interest to know what the integration time of the detector 
is. During the initial stages of evasion knowledge of the 
integration time is not critical. 

Graph l.D is a representative plot of the costate 
variables. The costate variables, in part, represent the 
importance of their corresponding state variable. Costate 
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variable X ^ corresponds to the x coordinate and X 2 to t ^ rie 
y coordinate. The graph shows that the evader, starting at 
state position (x = 60, y = 0) wants to increase y and 
values y slightly over x. That is, the evader values an 
increase in the lateral distance, y over an increase in 
time. This is true throughout the trajectory. 
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INFINITE INTEGRATION TIME 



RELATIVE MOVEMENT PLOT 
(WITH EVADER’S CONTROL VECTORS) 




For interpolation 
of control vectors: 
refer to graph 2. A. 2 
for velocity 
refer to graph 2. A. 3 
for heading 



VELOCITY SCALE-KNOTS 

-t- » f > -[ 



NOTE: 

1 . 

2 . 



GRAPH l.A 
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y coordinate -- miles 



ZERO INTEGRATION TIME 



relative movement plot 

O’/ITH EVADER'S CONTROL VECTORS ) 



SEARCH VELOCITY” 15.O KNOTS 




NOTE: For interpolation of 



VELOCITY SCALE— KNOTS 


1 # 


control variables: 
refer to graph 2.B.2 


for 


velocity 


| "t- ! 1 I | • r 1 | 

0 10 20 






2 . 


refer to graph 2.B.3 


for 


heading 



GRAPH l.B 
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Y COORDINATED — MILES 



INFINITE AND ZERO 

INTEGRATION TIME EVASION TRAJECTORIES 




GRAPH l.C 
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INFINITE INTEGRATION TIME 
COSTATE VARIABLES VERSUS X COORDINATE 



SEARCH VELOCITY =15.0 KNOTS 




X COORDINATES -- MILES 



GRAPH l.D 



Y COORDINATE — MILES 



INFINITE INTEGRATION TIME 



RELATIVE -MOVEMENT PLOT 



SEARCH VELOCITY =15.0 KNOTS 




GRAPH 2.A.1 
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VELOCITY IN KNOTS 



INFINITE INTEGRATION TIME 



15.0 

12.5 

10.0 

7.5 
5.0 

2.5 

0.0 



VELOCITY VERSUS X COORDINATE 



The curves are for the 
initial state positions 
shown in ( ) . 

SEARCH VELOCITY =15.0 KNOTS 




GRAPH 2. A. 2 
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PHI — IN RADIANS 



INFINITE INTEGRATION TIME 
COURSE HEADING VERSUS X COORDINATE 



The curves are for the 
initial state positions 
shown in ( ) . 




X COORDINATE -- MILES 



GRAPH 2. A. 3 
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NORMALIZED POWER DENSITY 
(unitless) 



INFINITE INTEGRATION TIME 

POWER DENSITY RECEIVED BY SEARCHER 
VERSUS X COORDINATE 

The curves are for initial 
state positions shown in 
( ). 



.008 



.006 
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000 



SEARCH VELOCITY =15.0 KNOTS 




X COORDINATE -- MILES 



GRAPH 2. A. 4 
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Y COORDINATES -- MILES 



ZERO INTEGRATION TIME 



RELATIVE MOVEMENT PLOT 



SEARCH VELOCITY - 15.0 KNOTS 




X COORDINATE — MILES 



GRAPH 2.B.1 



40 



VELOCITY IN KNOTS 



ZERO INTEGRATION TIME 



VELOCITY VERSUS X COORDINATE 



The curves are for the initial 
state positions shown in ( ) . 




X COORDINATE -- MILES 



GRAPH 2.B.2 
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PHI in RADIANS 



ZERO INTEGRATION TIME 



HEADING VERSUS X COORDINATES 



The curves are for the initial 
state positions shown in ( ) . 



SEARCH VELOCITY =15.0 KNOTS 




X COORDINATE -- MILES 



GRAPH 2.B.3 
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Y COORDINATES -- MILES 



INFINITE INTEGRATION TIME 
RELATIVE MOVEMENT PLOT 



SEARCH VELOCITY =10.0 KNOTS 




NOTE: Control variables may 

be obtained from 
Graphs 3. A. 2 and 3. A. 3 



GRAPH 3.A.1 
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VELOCITY -- KNOTS 



INFINITE INTEGRATION TIME 



EVASION VELOCITY VERSUS X COORDINATE 



The curves are for the Initial 
state positions shown in ( ) . 



Search velocity = 10.0 knots 




GRAPH 3 . A . 2 
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RADIANS 



INFINITE INTEGRATION TIME 



COURSE HEADING VERSUS X COORDINATE 



The curves are for the initial 
state positions shown in ( ) . 



SEARCH VELOCITY =10.0 KNOTS 




GRAPH 3. A. 3 
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NORMALIZED POWER DENSITY (unitless) 



INFINITE INTEGRATION TIME 

POWER DENSITY RECEIVED BY SEARCHER VERSUS X COORDINATE 



The curves are for the initial 
state positions shown in ( ) . 




GRAPH 3. A. 4 



60 

50 

40 

30 
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10 



ZERO INTEGRATION TIME 



RELATIVE MOVEMENT PLOT 




10 30 30 40 SO 60 

K COORDINATE — NILES 



NOTE: Control variables may be 

obtained from Graphs 3.B.2 
and 3.B.3. 



GRAPH 3.B.1 
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VELOCITY IN KNOTS 



ZERO INTEGRATION TIME 



VELOCITY VERSUS X COORDINATE 



The curves are for the initial 
state positions shown in ( ). 



SEARCH VELOCITY =10.0 KNOTS 




GRAPH 3.B.2 
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PHI IN RADIANS 



ZERO INTEGRATION TIME 



COURSE HEADING VERSUS X COORDINATE 



The curves are for the initial 
state positions shown in ( ), 



SEARCH VELOCITY =10.0 KNOTS 




X COORDINATE -- MILES 



GRAPH 3.B.3 
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APPENDIX A 



SOLUTION TO INFINITE INTEGRATION TIME PROBLEM 



The problem may be stated as. 



MAX J 




x 




dt 



subject to; 

dt = - V S + v sin 0 



= v cos 0 
x(t 0 ) = x 0 

y( t o ) = yo 

x(t X ) = x x 

y(ti) = y x 

where P(v) = 1. + 0.00005 : Velocity, v, and course 

heading, 0 , are functions of time. Velocity is measured 
in knots. The limits of integration should, in theory, be 
As stated tQ = 0 and infinity is approximated with a large 
number t-^. The symbols are as shown in Figure 1. 



dy 

dt 
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I 



V 



y 




The Hamiltonian function is now defined as 

H (x,y,v,0,t) = I (x,y,v,0,t) + A f (x,y,v,0,t) . 
In the case under discussion, 



H(x,y,v,0,t) = ~Z P + y y2 ’ + ll (' v s + v sin 0 ) +A 2 v cos 0 

( 1 ) 

The object is to maximize the Hamiltonian function at 
every point of time along the optimal trajectory by choice 
of the control variables, v and 0 . It is required that; 



o 

II 

X > 

CO 


Ai = 

1 d x 


B H = o 


A - - C> H 
2 3j 


AW - 0 . 


a 2 <h> - » 
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( 2 ) 



Now, from (1) and (2), 







2 x P(v) 
(x 2 + y 2 ) 2 



( 3 ) 




2 y P(v) 

(x 2 + y 2 )2 



(4) 



Also, 



d h 

3 v 


P'(v) 
x^ + y 2 


+ A^ sin 0 + A 2 cos 0 


(5) 


3 H 

30 " 


A x v cos 0 


- A 2 v sin 0 


(6) 



Setting equation (6) equal zero implies, 



^ 1 _ sin 0 

A 2 cos 0 

0 = tan -1 (A^AP * (?) 



From trigonometry and Figure 2 it can be seen that, 



and 



COS 0 



sin 0 



^2 

( A x 2 + A 2 2 ) % 

li 

( Ai 2 + A 2 2 



( 8 ) 

(9) 



52 




Figure 2 



Substitution of equations (8) and (9) into (5) results in; 
P'(v) _ A 2 + A 2 



* 2 + y 2 ( A1 2 + A 2 2 > 



*\/ A 1 2 + A; 



Therefore, P'(v) = (x 2 + y 2 ) V A i 2 +I 2 2 



( 10 ) 



Since P (v) = 1.0 + 0.00005 v , 

P'(v) = 0.0002 v 3 



(ID 



Equating (11) and (10) and solving for v results in 

V = (5000. (x 2 + y 2 ) *YAi 2 + A 2 2 ) 1/3 • (12) 

Equations (7) and (12) are the conditions for optimality, 
and must be satisfied at every point in time along the 
optimal trajectory. 

Since A^(t^) = = from (12) it follows that 

v(t^) = 0 . (13) 
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Since 



0 (tp = tan 



A.2<ti) 



is an indeterminate 



-1 



form a solution is sought in L'Hopital's rule. Rewriting (7), 



tan <p( t L ) = 






and applying L'Hopital's rule 



3 Ai^i) 

lim tan 0 (tp = — ^ 1— . (14) 

t->t x d A2^l) 

~dt 



Substituting (3) and (4) into (14) and canceling similar 
terms results in 



lim 

t-»t 



1 



tan 0 (tp 



_x 

y 



(15) 



Rewriting (15) the terminal condition for the control 
variable ft is 



0 (t0 = tan' 1 -2- 



(16) 



The objective is to find the values of the control variables 
at every point in time. Of principle interest are the 
initial velocity and heading. The solution now involves 
four simultaneous partial differential equations. The 
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initial conditions of variables x and y could easily be 
given, but only the terminal conditions are known for 
X ^ and ^ 2 * The solution is to reverse time and call 
the initial conditions terminal conditions and vice versa 
Reversing time changes the signs of the differential equa 
tions, but has no effect on the optimal conditions. 

The problem solution in backward time may now be 
stated as, 



o = 

X 


V g - v sin 0 


o 

y 


- v cos 0 


o 


2 x P(v) 


Ai " 


(X 2 + y 2 ) 2 


ii 

CM 


2 y P(v) 

( X 2 + y2)2 



where the "o" represents backward time. The initial 
conditions are; 

\ 

x(t ) = parameter 

\ state variables 

y(t ) = parameter 
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costate variables 



A^to) = o 

A 2 (t 0 ) = 0 



v(t Q ) = o 

0 (t Q ) = tan 1 (x/y) 



l control variables 
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APPENDIX B 



COMPUTER OUTPUT AND METHODS OF SOLUTION 

The computer routines given in this appendix were set 
to operate on the time sharing system (cp/cms) of an IBM 
360 computer. The graphical display package of a Tektronix 
4012 display terminal were used. Various termination state 
positions were attempted until an initial evasion range was 
selected. Since all the information for a particular 
trajectory was calculated at each execution, it was only 
necessary to change the call statements in the plotting 
package for different visual display. 
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FLOW CHART FOR INFINITE TIME SOLUTION 




T is in hours 
v is in knots 
x,y are in miles 
V g is searchers speed 
in knots 
P(v) = 1+ b v* 



NO 



IS Y LESS THAN 0 ? 



- YES 


STOP 




— 
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